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Abstract 

The present Chapter discusses methods by which topological Bloch bands 
can be prepared in cold-atom setups. Focusing on the case of Chern bands 
for two-dimensional systems, we describe how topological properties can be 
triggered by driving atomic gases, either by dressing internal levels with light 
or through time-periodic modulations. We illustrate these methods with 
concrete examples, and we discuss recent experiments where geometrical 
and topological band properties have been identified. 


1.1 Introduction 

Ultracold atoms constitute a promising physical platform for the preparation 
and exploration of novel states of matter EM. In particular, the engineer¬ 
ing of topological band structures with cold-atom systems, together with the 
capability of tuning interactions between the particles, opens an interesting 
route towards the realization of intriguing strongly-correlated states with 
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topological features, such as fractional topological insulators and quantum 
Hall liquids @j. 

This Chapter is dedicated to the preparation and the detection of topo¬ 
logical band structures characterized by non-zero Chern numbers [5]. Such 
Chern bands, which constitute the building blocks for realizing (fractional) 
Chern insulators mm , arise in 2D systems presenting time-reversal-symmetry 
(TRS) breaking effects. For instance, non-trivial Chern bands naturally ap¬ 
pear in the Harper-Hofstadter model |8j, a lattice penetrated by a uniform 
flux, where they generalize the (non-dispersive) Landau levels to the lat¬ 
tice framework. Additionally, Chern bands also appear in staggered flux 
configurations, such as in Haldane’s honeycomb-lattice model [9], or in lat¬ 
tice systems combining Rashba spin-orbit coupling and Zeeman (exchange) 
fields. 

The atoms being charge neutral, “magnetic” fluxes cannot be simply pro¬ 
duced by subjecting optical lattices to “real” magnetic fields. It is the aim of 
this Chapter to review several schemes that have been recently implemented 
in laboratories with the goal of realizing synthetic magnetic fields leading 
to Chern bands for cold atoms. The Chapter is structured as follows: Sec¬ 
tion 1.2 describes how the Chern number is related to physical observables 
defined in a lattice framework. In particular, it clarifies the link between re¬ 
cent Chern-number measurements performed in cold bosonic gases and the 
more conventional (electronic) quantum Hall effect. Section 1.3 reviews a 


few lattice models, and it relates their properties to the well-known Landau 
levels of the continuum. The main Section 1.4 describes diverse schemes by 
which Chern bands can be prepared and probed for cold atoms; a special 
emphasis is placed upon experimentally realized schemes. The last Section 
11.51 is devoted to final remarks and discussions. 


1.2 Bloch waves, the anomalous velocity 
and the Chern number 

The Chern number - a topological invariant v c h G Z classifying fibre bun¬ 
dles - naturally enters the description of particles moving in two-dimensional 
lattices, where it offers an elegant interpretation for the (anomalous) quan¬ 
tum Hall effect mmm- This Section relates this topological invariant to 
physical observables and discusses methods to measure it in experiments. 
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1.2.1 Bloch bands, the anomalous velocity and the Chern number 


The general problem of a particle subjected to a 2D space-periodic potential 
U(r+a) = U (r) starts by invoking Bloch’s theorem, which stipulates that the 
eigenstates of the system can be decomposed as Vt\fc( r ) = exp(ir • k) u Xk (r ), 
where u Xk (r) has the periodicity of the lattice u Xk (r + a) = u Xk {r) and 
where k = (k x ,k y ) is the quasi-momentum. The associated eigenenergies 
E x (k) display bands, labeled by the index A, over the first Brillouin zone 
(FBZ) of the quasi-momentum k. In the absence of additional potentials, a 
state u\k in a given Bloch band A with quasi-momentum k is characterized 
by the averaged velocity v x (k) = (u Xk \v\u X k) = (1 /h)dE x (k)/dk, where v 
denotes the velocity operator. 

Subjecting the lattice system to a constant force F generates Bloch oscil¬ 
lations. Indeed, restricting the dynamics to ID for now, and considering for 
simplicity a semi-classical approach, the system is described by the equations 
of motion 

hx c (t) = Hv\(k c ) = > hk c (t) = F, (1.1) 


where x c [t) and hk c {t) denote the center-of-mass position and momentum 
of a wave packet prepared in the lowest band E\ (i.e. A = 1), see Fig. 0 
The equations of motion (1.1) are valid when the force F is weak enough to 
preclude any inter-band transitions. 


t = 0 t 


E^k) 




k° c +Ft/h k 



Figure 1.1 Bloch oscillations for a wave packet initially prepared in the 
lowest band E\, with center-of-mass position x® and momentum fc °. T he 
trajectory x c (t) results from the semi-classical equations of motion (1.1). 


When considering 2D lattices, an interesting effect adds to the standard 
Bloch oscillations, which involves the Berry curvature of the band |12j. In¬ 
deed, applying a force along a given direction, e.g. F = F y l y , modifies the 
averaged velocity v\(k) along the transverse (x) direction according to 


v x x(k) 

n x (k) 


dE x (k ) 
hdk x 


F,„ 


n x (k), 


i {(dk x u X k\d ky u Xk ) - (d ky u Xk \d kx u Xk )) , 


( 1 . 2 ) 

(1.3) 
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where Q\(k) denotes the Berry curvature of the band A. Thus, the velocity 
in a state u\(k) has two contributions: the usual band velocity, responsible 
for Bloch oscillations [Fig. O’ and the so-called anomalous velocity, which 
can produce a net drift transverse to the applied force. As for Eq. (1.1), the 


general result (1.3) assumes the absence of inter-band transitions (i.e. weak- 
force regime). 


The anomalous (transverse) velocity in Eq. (1.3) can be isolated and ob¬ 


served experimentally by canceling any contribution from the band velocity, 
which can otherwise dominate. This can be achieved by comparing trajec¬ 
tories for opposite forces, ±E y , noting that the anomalous velocity changes 
sign under reversal of the force m- Alternatively, the average anomalous 
velocity can be isolated by uniformly populating the bands, namely by av¬ 
eraging the velocity in Eq. (1.3) over the entire FBZ, since 


/ FBZ 


(dE\(k)/dk x ) d 2 k = 0, 


(1.4) 


due to the periodicity of the energies in k- space. 

We now compute the total transverse velocity in the case of uniformly 


populated bands using Eq. (1.3); explicit physical implementations will be 


discussed in the next Sections 1.2.2|[LT3 In the following, we consider a 
general square lattice system of size A sys t = L x xL y , characterized by a unit 
cell of size A ce n = d x a x d y a, where a denotes the primitive lattice spacing 1 . 
The number of states within each band is fV s tates = A sys t/A ce n. We write the 
total number of particles as Wot = where is the number of 

particles occupying a given band A. We now make the assumption that each 
band is populated homogeneously, so that the average number of particles 
in a state u\k is uniform over the FBZ, and it is given by 


P W (k) = p™ = N^/N s tates. 


(1.5) 


Using Eqs. (1.3) and B we obtain the total averaged velocity along the 
direction transverse to the force 


Wot =xy A) XM( fc ) 

A k 

1 

2tt 


F yA ce U^AtA)^), 


h 


A 


u ch = 7W «A(fc) Ak x Ak y 


7^ / n x (k)d 2 ke z, 

27T J FBZ 


( 1 . 6 ) 

(1.7) 


where A k XtV = 2-k/L x ^. The latter equations reveal that the total velocity 

1 For instance, d x =d y /2 = 1 for the brick-wall lattice m, and d x =d y = 2 for the magnetic unit 
cell of a square lattice penetrated by a uniform flux <I> = 7r/2 per plaquette m- 
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is related to the quantities which converge towards the Chern num¬ 
ber of the bands A when taking the thermodynamic limit L x>y —> oo. The 


Chern number is an integer, obtained by averaging the Berry curvature 
n x {k) over the FBZ [Eq. (1.7)]; it is a topological invariant, meaning that 
remains a constant as long as the spectral gaps to other bands do not 
vanish, see Refs. mm • As announced above, any contribution from the 
band velocity cancels under the homogeneous-population condition (1.5), as 
a direct consequence of Eq. (1.4). 


The following Sections 1.2.2 1.2.3 discuss two different physical realiza¬ 


tions that revealed the Chern numbers v. 


(A) 

ch 


homogeneous population of energy bands [Fig. 1.2 


in experiments, through the 
These Sections aim to 


clarify the link between the quantum Hall effect 10. ITT! [5], as observed 
in electronic systems since the 1980’s, and the Chern-number measurement 
recently performed with ultracold bosonic atoms M- 


• Filled band of fermions 


• Thermal gas 




e 2 



Ei ^ Rband ksT 


P — Aiol /Abates 1 


Figure 1.2 Two realizations of uniformly populated energy bands, suitable 
to reveal the Chern number in experiments, (left) Considering fermions at 
zero temperature, the band is perfectly filled by setting the Fermi energy 
within the spectral gap. (right) A system of bosons uniformly populate the 
band when the temperature is large compared to the bandwidth Wband, 
but small compared to the gap. 


1.2.2 Fermions, the quantum Hall effect and the TKNN formula 

The first situation that we consider is a 2D non-interacting polarized Fermi 


gas at zero temperature [Fig. 1.2 . Setting the Fermi energy Ep within a 


spectral gap naturally leads to a perfect filling of the bands E\ < Ep located 
below the gap: the average number of particles in a state u x (k), Eq. (1.5), 
is exactly /N st & tes = 1 for E\ < Ep. Setting the latter condition 
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into Eq. (1.6) yields 


Vtot = “ 


FyA s y s t 

h 


E 

E\<Ey 


/A) 

ch > 


( 1 . 8 ) 


which indicates that the total velocity of the Fermi gas is directly related to 
the sum of Chern numbers associated with populated bands. 

In systems presenting time-reversal symmetry, the Berry curvature satis¬ 
fies £l\(k) = in which case all the Chern numbers = 0, see 

Eq. (1.7). Hence, observing the transverse drift associated with v x ot requires 
a system without time-reversal symmetry, which is the case, e.g., in elec¬ 
tronic systems subjected to magnetic fields [5]. In these electronic setups, 


the transverse transport predicted by Eq. (1.8) is measured through the Hall 


conductivity a xy , which relates the electric field E y to the current density 
j x = ev* ot /A s y S t, where e is the electron charge. Using these definitions to¬ 


gether with Eq. (1.8), we recover the well-known TKNN formula [10] for the 


electric Hall conductivity 


jx — CT X yE y , 


°y x ~ h ^ 

E\<E f 


'A 


= — a 


xy, 


(1.9) 


where we introduced the electric field E y = F y /e acting on the electrons. 


The TKNN formula (1.9) expresses the fact that transport measurements 
directly reveal the Chern numbers in electronic systems, through the 
quantization of the Hall conductivity. In particular, the quantized value 
a yx = (e 2 /h)x (integer) remains constant as long as the Fermi energy stays 
in an open gap mm- Note that populating a Chern band with a Fermi 
gas is not limited to electronic systems, as it could also be performed by 
trapping fermionic atoms in an optical lattice subjected to artificial magnetic 
fields [16]. 


1.2.3 Thermal Bose gas and the center-of-mass displacement 

Let us now consider a radically different configuration: a thermal gas of non¬ 
interacting polarized bosons, whose temperature is large compared to the 
bandwidth Uband of the lowest band E\ , but small compared to the spectral 
gap A above it [Fig. |1.2| . In this case, the average number of bosons in 
a state of the lowest band is homogeneous but density-dependent, pA) = 
= ^Vtot/-^states / 1) an d P ^ = 0 for A > 1. The homogeneity of the band 
filling, which is a reasonable assumption when the lowest band presents a 
large flatness ratio A/ITband ^ 1, can be tested in cold-atom experiments 
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using band-mapping techniques mug. Setting the particle filling condition 


into Eq. (1.6) yields 


«tot = - 


FyA ; 


h 


■syst (1) 

P v ch> 


P — -Wot /-^states 1 1 


( 1 . 10 ) 


which differs from the Fermi-gas result in Eq. (1.8) in that only the lowest 


band contributes, but also through the additional density dependence. 

The Bose gases considered in cold-atom experiments are charge neutral. 


However, in analogy with the quantum Hall effect discussed in Section 1.2.2 


one could consider a transport measurement relating the particle current 


density j x = uf ot /H syst to the applied force F y . In this case, using Eq. (1.10), 
we find that the analogue of the Hall conductivity would read 


jx — ®xyFyi &xy — (1 P ^ch • 


( 1 . 11 ) 


In contrast with the TKNN formula ( |1.9| ), the relation between the mea¬ 
sured transport coefficient a xy and the topological Chern number involves a 
density-dependent factor p/h. Thus, in such an experiment, identifying the 
Chern number of the lowest band would require to combine simultaneous 
transport and density measurements, which constitutes a severe drawback of 
this approach. To overcome this difficulty, one can consider another physical 
observable: the center-of-mass (CM) displacement of the gas. Indeed, using 


Eq. (1.10), we find that the transverse velocity of the CM is given by 


V CM ~ v tot/Ntot — — 


Py-^cell 

h 




( 1 . 12 ) 


Since both the unit cell area ^4 C eii and the strength of the applied force F y 
can be determined with precision, the transverse displacement of the CM 
A®CM(t) = VQ M t offers a direct measure of the Chern number of the low¬ 
est band. This Chern-number-measurement was successfully implemented in 


Munich m in 2014, see Section [L4| 


1.3 Models: Harper-Hofstadter bands vs. Haldane-like bands 

We now turn to the presentation of two lattice models that have been re¬ 
cently implemented with cold atoms setups and that lead to a topologically 
non-trivial lowest band. For the sake of completeness, we start our discussion 
with a short reminder of the well-known Landau levels, which characterize 
the quantum motion of a free particle in a uniform magnetic field. This 
Landau level structure, with its non trivial bulk topological features and the 










associated edge currents, is crucial to explain Quantum-Hall-type phenom¬ 
ena, both in the integer and the fractional cases. 


1.3.1 Landau levels 

The Landau level spectrum m emerges when one looks for the eigenstates 
of the Hamiltonian H = (p — qA{r)) 2 /2m, where p (resp. r) denotes the 
canonical momentum (resp. position) operator for a particle of charge q and 
mass m. The vector potential A(r), defined up to a gauge transformation, 
satisfies VxA = B, where B is uniform. In order to derive the spectrum of 
H , we note that it can be written as H = (n^+n^)/2m, where we introduced 
the kinetic momentum II = p — qA(r). The kinetic momentum operators 
satisfy the simple commutation relation [H x , H y ] = ihqB. The corresponding 
operator algebra is thus formally equivalent to that of a harmonic oscillator 
where H = (X 2 + P 2 )/2 with [A', P] = i. We infer from this equivalence that 
the energy levels form an equidistant set huj c (n + 1/2), where u c = qB/m 
denotes the cyclotron frequency and n is a nonnegative integer. Each Landau 
level has a macroscopic degeneracy S/(2tt£ 2 ), where S is the area of the 
sample and £ = (h/qB) 1 ^ 2 is the magnetic length. One thus obtains a band- 
like spectrum as for a particle in a periodic lattice, each band being infinitely 
narrow. Each Landau level has a Chern number equal to 1, so that the 
notion of anomalous velocity is also relevant here. Actually it has in this 
case a simple classical interpretation in terms of Hall current. It is indeed 
well known that when a particle in a magnetic field along z is acted upon 
by an additional uniform force F y , its motion consists in the combination 
of the circular cyclotron motion and a uniform translation motion at the 
constant velocity v x = F y /(qB)\ the latter is nothing but the anomalous 
velocity described above. 


1.3.2 Lattices and tight-binding models 

We now turn to the Harper-Hofstadter problem [18L IH], which is the trans¬ 
position of the Landau problem to a discretized space. It models the motion 
of a charged particle in a square two-dimensional lattice of period a, normal 
to the uniform magnetic field (Fig |1.3| (a)). This motion is described using 
the single band approximation, assuming that only nearest neighbour cou¬ 
plings are relevant (Hubbard model). The quantum state corresponding to 
the particle localized on the lattice site r = a(m,n) is denoted |r), with 
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Figure 1.3 Two lattice models leading to topologically non-trivial bands, 
(a) The Harper-Hofstadter model: a square lattice operating in the tight- 
binding regime is placed in a uniform field, with the flux $ across each unit 
cell of the lattice [TBj, [5]. (b) The Haldane model: a periodic honeycomb 
lattice with nearest-neighbour couplings A -f-)- B (continuous lines) corre¬ 
sponding to a real tunneling matrix element, and next-to-nearest neighbour 
couplings A o A and B o B (dotted lines) corresponding to a non-real 
tunneling matrix element. The arrow directions encode the sign of the ar¬ 
gument of this matrix element [9]. 


m, n integers. The Hamiltonian reads in zero magnetic field 

Ho = ~J ^2 (!-13) 

{ r,r') 


where a r annihilates a particle in the state |r) and J > 0 is the tunneling 
matrix element. The eigenstates of this Hamiltonian are the Bloch func¬ 
tions |V>fc) = Yl r elk r \ r ) w ith energy E(k) = —2 J [cos (k x a) + cos (k y a)], 
corresponding to an allowed band of width 8 J. 

The presence of the magnetic field is taken into account in the single¬ 
band approximation by assigning a complex value to the tunneling matrix 
elements between adjacent sites. The individual phases of these matrix el¬ 
ements have to be chosen such that the total Aharonov-Bohm phase m 
accumulated along a closed contour is equal to 27 t< 1>/<f>o, where is the flux 
of the magnetic field through the contour and <l>o = h/q is the flux quantum 
BM- Because of gauge invariance, there are an infinite number of choices 
for individual phases matching this prescription. Here we choose the Lan¬ 
dau gauge, which amounts to taking zero phase (a real matrix element) for 
tunneling along the y direction, and the y-dependent phase 27 r cm for the 
tunneling along the x direction, where we introduced the dimensionless flux 
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parameter a = ^/^o- The Hamiltonian is now 

H a = — (®L,n+l“™,n + e ^ln+l,n^rn,nj + H.C. (1-14) 

m,n 

which is periodic in a with period 1. The term exp(i27rcm), which captures 
the effect of the magnetic field, is generally referred to as the Peierls phase- 
factor mm- 

The single-particle spectrum plotted as a function of a acquires a fractal 
structure known as the Hofstadter butterfly. The calculation of this spectrum 
for an arbitrary value of a is a difficult mathematical task. If one restricts to 
rational values a = p/q (jp and q relatively prime integers), the problem is 
simpler since one recovers a two-dimensional periodic problem, with a period 
a along x and qa along y. The corresponding unit cell (the so-called magnetic 
cell) now contains q sites, and the energy band of width 8 J for $ = 0 split in 
q non-overlapping subbands (two adjacent subbands may touch each other 
via Dirac points). In general each subband has a nonzero Chern index, which 
can be obtained via the solution of a Diophantine equation m, whereas the 
sum of all Chern indices over the q subbands is null. In the particular case 
a = 1/q, the Chern index of the lowest band is 1, so that this band is 
topologically equivalent to the lowest Landau level found in the absence of 
a lattice. 

In the Harper-Hofstadter problem, the presence of the uniform magnetic 
field breaks the translational symmetry of the initial lattice. Haldane pro¬ 
posed in 1988 another model [9], in which a non-trivial band topology could 
appear without any modification of the lattice unit cell nor breaking of the 
translational symmetry. The starting point is a graphene-like honeycomb 
lattice with nearest-neighbour couplings (Fig. |1.3| (b)). The unit cell thus 
consists of two equivalent sites, denoted A and B, with tunneling from sites 
A (resp. B) to the three neighbouring sites B (resp. A) with equal amplitude. 
At this stage the single-particle spectrum consists of two subbands, touch¬ 
ing each other in Dirac points. Haldane’s crucial insight was to add non-real 
next-to-nearest-neighbour (NNN) couplings, i.e. A —> A and B —>• B, which 
break time-reversal symmetry. By contrast to the Harper-Hofstadter model, 
these couplings are constant over the whole lattice, and could in principle 
be created by a staggered magnetic field with a zero-flux through the hon¬ 
eycomb unit cell. These additional couplings lift the initial degeneracy at 
the Dirac points, and the two subbands are now separated by a gap with 
non-zero opposite Chern indices, +1 and — 1. With the lowest band filled 
with spinless non-interacting fermions, Haldane’s model constitutes a proto¬ 
type of a Chern insulator that goes beyond quantum Hall setups. It played 
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a key role for the subsequent discovery of time-reversal invariant topological 
insulators, where spin-orbit coupling replaces the complex next-to-nearest- 
neighbour tunnel coupling. The optical flux lattice concept that will be pre¬ 
sented latter in this Chapter is also directly related to Haldane’s model, 
since it provides topologically non-trivial energy bands for configurations 
where the atom-light interaction is periodic. 


1.4 Implementation: driving atoms into topological matter 


This Section describes several schemes realizing Chern bands in cold-atom 
systems. For the sake of clarity, this Section makes the distinction between 
methods based on a tight-binding approach [Section |1.4.1 and those ap¬ 
plicable in the continuum or weak-lattice regimes [Section 1.4.2 . A special 


emphasis is set upon recent theoretical works and experimental implemen¬ 
tations, which led to unambiguous signatures of topological properties asso¬ 
ciated with non-zero Chern numbers in 2D optical lattices. 


1.4-1 Schemes based on a tight-binding approach 

Non-interacting cold atoms moving in a deep optical lattice [I] are well de¬ 
scribed by the single-band tight-binding Hamiltonian in Eq. ( 1.13| ). Starting 
with this topologically-trivial tight-binding band, it appears that a nat¬ 
ural way to generate Chern bands consists in implementing the Harper- 
Hofstadter or the Haldane model [Section 1.3.2 , which can be achieved by 


controlling the tunneling matrix elements within the optical-lattice setup, 
see Eq. (1.14). Indeed, as was discussed in Section 1.3.2, Peierls phase fac¬ 


tors can be associated with (local) magnetic fluxes, which break TRS and 
potentially lead to Chern bands. We now review several methods by which 
Peierls phase factors and magnetic fluxes can be engineered using different 
aspects of cold-atom technology, see Fig. |1.4[ 


Using internal states of the atoms 

A natural way to induce complex tunneling matrix elements in optical lat¬ 
tices is to exploit the spatial-dependence of the optical phase in photon- 
assisted tunneling mmmm- Consider two internal states of an atom, 
denoted \g) and |e), respectively trapped in two independent (state-dependent) 
optical lattices V 9ie along the x direction, see Fig. |1.4[ b). For large enough 
lattice potentials, and in the absence of coupling between the two states, 
the hopping is completely inhibited. Adding a coherent coupling, with wave 
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vector q and frequency uj ge matching the energy difference between the in¬ 
ternal states, effectively activates tunneling processes between the two sub¬ 
lattices. The tunneling matrix elements between two nearest-neighbouring 
lattice sites r g and r e are of the form J(r g ) = l 7 e frexp(iq • r g ), where the 
amplitude is proportional to the coupling’s Rabi frequency and to the 
overlap between Wannier functions defined at the two sites involved in the 
process mmmm- Hence, space-dependent Peierls phase-factors, and the 
corresponding magnetic fluxes penetrating lattice plaquettes, can be con¬ 
trolled by tuning the coupling’s wave vector q. This method can be applied 
to square optical lattices, leading to an effective Harper-Hofstadter Hamil¬ 
tonian (1.14), but also to triangular/honeycomb lattices in view of realizing 
the Haldane model mmm- The scheme is general and can be applied 
to diverse atom-light configurations: for instance, it may involve two hy- 
perfine states in the ground-state manifold coupled through a (two-photon) 
Raman coupling [20], or a one-photon coupling between a ground state and 
a long-lived excited state [ 21 ] • This method has not yet been implemented 
experimentally. 

The scheme described above concerns the coupling between atoms living 
on a 2D optical lattice. However, it is intimately related to the concept of 
“synthetic magnetic fields in synthetic dimensions” E3, which has been im¬ 
plemented at LENS [28j and NIST [29] in 2015. Here, atoms in M internal 
states live on a standard ID optical lattice, characterized by real-valued tun¬ 
neling matrix elements. Atom-light coupling then drives on-site transitions 
between the M internal states, such that each transition m —> m + 1 is ac¬ 
companied with a space-dependent phase q ■ r (i.e. a momentum transfer). 
Interpreting these internal-states transitions as hopping processes along a 
synthetic (internal-state) dimension, which is spanned by m = 1 ,...,M, 
this system is found to be equivalent to the effective Harper-Hofstadter sys¬ 
tem described above. The propagation of chiral edge-states, an unambiguous 
signature of synthetic magnetic fluxes in 2D lattices, has been experimen¬ 
tally measured in these experiments 


Superlattices and resonant time-modulations 
Complex tunneling matrix elements can also be engineered without consid¬ 
ering the internal structures of the atoms. In direct analogy with the scheme 
described above, this can be achieved by imposing an energy offset between 
neighboring sites of an optical lattice and modulating the system with a 
resonant time-periodic modulation of the on-site energies, see Fig. |1.4| (d). 
In this case, the Peierls phase-factors are directly related to the phase of 
the modulation. To be explicit, let us consider a minimal lattice system: two 



Preparing and probing Chern bands with cold atoms 


13 


(a) Coal: induce Peieris phase factors (b) Using internal states of the atoms 



Figure 1.4 (a) Effective magnetic fluxes can be created in deep optical 
lattices by inducing complex tunneling matrix elements, (b) Method us¬ 
ing state-dependent optical potentials V g ^ e , trapping atoms in two internal 
states | < 7 ) and |e), combined with a resonant atom-light coupling with fre¬ 
quency uj ge = (E e — E g )/h. (c) Shaking a lattice, with an arbitrarily large 
driving frequency ui. (d) Method using a superlattice with energy offset 
A, large compared to the bare hopping energy, combined with a resonant 
onsite energy modulation with frequency ui = A/h. 


lattice sites (a and 6 ), separated in energy by a large offset A, and subjected 
to a resonant modulation acting on the lower site (a) only. We write the 
tight-binding Hamiltonian in the form 

H(t) = —J + tfa'j + A b^b + k cos(tuf + fyaJa, (1.15) 

where J is the bare hopping amplitude between the sites. In the large- 
frequency regime A = Huj 3> J, the long-time dynamics resulting from 
the time-dependent Hamiltonian H(t) is found to be captured by a time- 
independent (effective) Hamiltonian, which in this case, takes the simple 
form EDI ED EH S3 El] 

H e ff = — J exp (i(f>) a)b + h.c., (1.16) 


where J\ denotes the Bessel function of the first kind. According to Eq. 
(1.16), the tunneling between the sites is effectively restored and the cor¬ 
responding tunneling matrix elements include a Peieris phase-factor related 
to the phase of the modulation </>, which can be made space-dependent. The 
result in Eq. (1.16) constitutes the building blocks for the generation of mag¬ 
netic fluxes in time-modulated optical lattices, as implemented in Munich 
[35l S3 US] and at MIT [37]. The following of the Section describes these 
schemes in more detail. 

In order to treat time-modulated optical lattices, we first introduce a set of 
useful equations, which offer a powerful theoretical framework to analyse the 
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physics of time-dependent problems. Let us consider a general Hamiltonian 
of the form H(t) = Ho + V(t), where V(t + T) = V(t) and where T = 2ir/u 
is the period of the driving. In the following, the period T is considered to 
be small compared to all characteristic time-scales [38]. We are interested 
in describing the dynamics of a initial state ip (to) under the driving, which 
is assumed to start at time to. Formally, the state at a given time ip(t) is 
obtained by acting on the initial state ip (to) with the time-evolution operator 
U(t\ to), which can be partitioned as 139, 38] 

U(t;t 0 ) = exp [~ik(t)} exp[-*(f - t 0 )H e s/K } exp[LR'(f 0 )], (1-17) 


where H c g is a time-independent (effective) Hamiltonian describing the long¬ 
time dynamics, and where the “kick” operator K(t) captures the micro¬ 
motion and the effects related to the initial phase of the modulation. The 
dynamics is thus completely captured by the operators H e g and K(t), which, 
according to Refs. [39, |3SJ !33], can be calculated systematically through 
a perturbative expansion in powers of 1/u. Considering a single-harmonic 
modulation, V(t) = V exp (iut) + h.c., these operators are approximatively 
given by |38] 

H eS = H 0 + (1 /hu)[V,V*] + 1 /(V2Hu) 2 , Ho],V^] + h.c.) + 0( 1/a; 3 ), 

K(t) = (l/ifrw) (Vexp(iuit) — h.c.^ + 0(1/J 2 ). (1.18) 


In some cases, the infinite series in Eq. (1.18) can be (partially or totally) 
resummed; this allows convergence of the expansion even in the strong¬ 
driving regime [38l IM1 [4Cl]. The equations (1.17)-(1.18) offer a systematic 
way to build the time-evolution operator, which can be directly applied to 
arbitrarily complicated systems. 

Having defined a theoretical framework to treat time-dependent problems, 
we now apply it to time-modulated optical superlattices. The Hamiltonian 
is taken in the form H(t) = Ho + V(t), where the static part consists of a 2D 
tight-binding Hamiltonian with an additional superlattice potential directed 
along the x direction 


Ho — + h.C.^ + A s(m)a\ nn a rrlin , 

m,n m,n 

(1.19) 

where a m ^ n creates an atom at lattice site r = a(m,n), a is the spacing 
and (m, n) are integers. The superlattice function s(m) is assumed to create 
energy offsets A J x between all neighboring sites, i.e. s(m+l) — s(m) = ±1 
for all m. The tunneling is then restored by applying an on-site time-periodic 






Preparing and probing Chern bands with cold atoms 


15 


modulation of the form 

V(t) = V exp(iu;t) + h.c., V = n Y, v(m, n a m , n - (1.20) 


m,n 


Here, the complex numbers v{m 1 n ) capture the space-dependence of the 
modulation, which turns out to be crucial, and we impose the resonance con¬ 
dition uj = A/h. The long-time dynamics of the modulated 2D lattice [Eqs. 


(1.19)-(1.20)] is captured by the effective Hamiltonian H e g- in Eq. (1.18), 
which after partial resummation of the series yields a Harper-Hofstadter- 
like Hamiltonian 


h oS -Yj x ( 


m n)p ±ie( ' m ’ n '>h^ h 

m ’ a ) e a m+l,n a m ,i 


+ l 7 2/ (m,n)aJ nr)+1 a min +h.c., (1.21) 


where 9(m,n ) = arg[u(m + l,n) — v(m,n)\, and where the sign ± depends 
on whether the hopping m —> m + 1 starts from a low- or a high-energy site 
of the superlattice potential s(m). Hence, the Peierls phase-factors in the 
effective Hamiltonian (and the corresponding magnetic fluxes per plaquette) 


depend on the time-modulation V(t), through v(m,n) [Eq. (1.20)], but also 


on the static superlattice potential s{m ) [Eq. (1.19)]. The effective hopping 
amplitudes J XiV are a priori space-dependent (see [Tl] and the discussion in 
my, they are explicitly given by 

J x (m, n) = J x Ji[K 0 \S x v(m, n)|], J y {m, n) = J y J 0 [K 0 \S y v(m, n)|], 

where I\q =2 k/Hco and 5 x _ y v denote finite-difference operations along the x 
and y directions, e.g. S x v(m,n) = v(m + 1 ,n) — v(m,n). The kick operator 
can be obtained along the same line; its impact on physical observables (e.g. 
momentum distributions) was discussed in Ref. 


Based on the result (1.21), two different features of the system should 


be designed so as to eventually realize the Harper-Hofstadter Hamiltonian 
( |1. 14 ): the superlattice function s(m) and the time-modulation function 
v(m,n). The main difficulty to achieve the uniform-magnetic flux config¬ 
uration is to handle the sign of the phase ±0(m, n), which can typically lead 
to staggered flux patterns in arbitrary superlattices |35j- To solve this issue, 
the Munich [36j and MIT [37] teams first combined a Wannier-Stark ladder 
potential generated by a magnetic-field gradient, i.e. s(m) = m, with a sim¬ 
ple moving potential, v(m, n ) = exp (iq ■ r), created by a single pair of laser 
beams. In this configuration, a uniform flux 27ra = q y a, together with ho¬ 
mogeneous hopping amplitudes J x ,y, were achieved. By setting q y = n/2a, 
the Munich team obtained a topological band structure, whose lowest band 
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had a Chern number v i = 1. However, 


able for the Chern-number measurement described in Section 1.2.3 


this setup was found to be unsuit- 

due to 


instabilities arising from the linear gradient s(m) = m. In order to achieve 
uniform flux with an all-optical superlattice [MIES], the Munich team then 
developed a novel setup, based on a two-site superlattice, s(m) = (—l) m , 
combined with a more sophisticated time-modulation of the lattice 3 . The 
latter was induced by two pairs of laser beams, so as to generate a uniform 
flux over the entire lattice pa eh, hence providing a stable platform to 
measure the Chern number, as we now describe. 

The Chern-number measurement m was achieved by loading bosonic 
atoms into the lowest band of the Harper-Hofstadter spectrum for a syn¬ 
thetic magnetic flux a = 1/4. The transverse drift of the cloud [Eq. (1.12)] 
was then detected in-situ , as a response to a weak optical gradient applied 
along the y direction. For short times, the transverse motion of the cloud 
was found to be linear, in agreement with the prediction in Eq. (1.12). For 
longer-times, heating processes, which were found to be independent of the 
applied force, promoted atoms to higher bands and lead to a saturation 
of the cloud’s transverse drift. A careful analysis taking into account the 
dynamical repopulations of the bands revealed the expected Chern number 
with a precision at the 1 % level. This Chern-number measurement was a di¬ 
rect probe for the topological order associated with the effective bulk energy 
bands, which resulted from the modulated optical superlattice. 


Off-resonant shaken optical lattices 


Before ending this Section on tight-binding realizations of Chern bands, let 
us briefly describe another promising strategy for creating and observing 
topological properties in optical-lattice systems. Similarly to the schemes 
discussed above, this method is based on time-periodic driving of the lattice, 
however, the driving frequency is now considered to be off-resonant with 
respect to any energy separations in the problem, see Fig. |1.4| )c). Let us 
compare this idea with the resonant-modulation approach, by re-writing 


the two-site minimal model in Eq. (1.15) as 


H(t) = — J ( 'd^b + b^a'j + «cos(u;t + (f)a)a, (1-22) 


and by computing the corresponding effective Hamiltonian (1.18) EH EH] 
H e ff = — J Ja{n/hw) (a^b + Pa'j . (1.23) 


3 Combining a two-site superlattice to a simple moving potential with v(m, n) = exp (iq ■ r ) 
leads to a staggered flux configuration, which is associated with zero Chern numbers |35l . The 
time-modulation used in El allowed to rectify the flux, by individually addressing successive 
links with two independent pairs of lasers. 
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Contrary to the resonant-modulation case [Eq. (1.16)], the driving phase 
(j) does not contribute to any Peierls phase-factor. In fact, a non-sinusoidal 
driving is required to generate effective complex tunneling matrix elements. 
This was demonstrated in Hamburg in 2012, where a Peierls phase factor was 
shown to appear in a shaken ID optical lattice (45] . The shaking strategy 
was also applied in 2D triangular lattices by the same team, in view of 
studying frustrated magnetism in optical lattices 116,3T]. While these setups 
indeed produce local magnetic fluxes, their corresponding band structure is 
associated with zero Chern numbers. 

In 2014, Jotzu et al. achieved the shaking of a 2D honeycomb lattice in a 
circular manner [38]. Similarly to the case of rotating traps [Section |1.4. 2 


this fast circular motion induced a chirality within the system, which for¬ 
mally breaks time-reversal symmetry. In fact, the effective Hamiltonian was 


shown to be equivalent to the Haldane model presented in Section 1.3.2 


the circular shaking effectively induces complex NNN tunneling matrix el¬ 
ements. This opens a bulk gap in the honeycomb-lattice (Dirac) spectrum, 
and it generates Chern bands. The Berry curvature of the bands was probed 
in the system |48j, through the observation of an anomalous velocity. 


1.4-2 Schemes in the continuum 

A natural route to the formation of topological bands with non-zero Chern 
number is to simulate the orbital effect of a uniform magnetic field. The 
atom should move continuously through space, subject only to this uniform 
field. This leads to the Landau level spectrum: highly degenerate bands, 
each with unit Chern number. There exist several ways to achieve this goal 
for cold atoms, at least to a good approximation. (Deviations can arise from 
non-uniformity of the field and/or the presence of additional potentials.) 

Rotation 

One very direct approach is to cause the atomic cloud to rotate, for example 
by applying a rotating deformation and bringing the system to equilibrium in 
a frame of reference rotating at angular frequency ft [U 09]. In this rotating 
frame an atom of mass m moving with velocity v experiences a Coriolis force 
F c = 2 mv x ft. This plays the same role as the Lorentz force on a charged 
particle in a uniform magnetic field, = qv x B. Equating coefficients, 
qB = 2?nfi, one deduces an effective flux density = q\B\/h = 2m\ft\/h. 
For a Bose-Einstein condensate, this manifests itself as a lattice of quantized 
vortex lines of areal density n^. Such vortex lattices have been observed, 
and their properties studied in detail, in experiments dating back over 15 
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years E23EDC2 [53]. For a circularly symmetric harmonic confinement of 
frequency ujq, mechanical stability under the centrifugal force sets an upper 
limit on rotation rate |fi| < ujq. At the point of balance, |fi| = ujq, the 
spectrum for the 2D motion is exactly that of Landau levels. This condition 
|n | ~ ujq has been reached experimentally to within 1% [54] . This rotation 
technique can also be implemented with optical lattices, see e.g. [55ll56l 5?]. 


Raman dressing 

An important feature of the physics of cold gases is that atoms (or molecules) 
can be prepared in different internal states (spin states or electronic excited 
states) which can have long lifetimes. It is then possible to use optical fields 
to drive atoms into coherent superpositions of these internal states. Opti¬ 
cal dressing of internal degrees of freedom provides very powerful ways to 
generate artificial gauge fields and to form optical lattices with topological 
bands. 

The physics underlying the use of internal states to generate artificial 
gauge fields can be understood in terms of a Berry connection in real space 
[58 , 59j. Consider an atom with two internal states, a = 1,2. We denote 
the eigenstates at position r for these internal states by |l) r and |2) r . In 
the presence of (optical) fields that couple the internal states, and in the 
rotating wave approximation, the local energy eigenstates are dressed states 
| n) r = ®n,i(f)|l)r + a n ,2 {f )|2)r• The coefficients a n ^ a (r) are determined by 
the fields at position r, so vary in space. Hence, the adiabatic motion of 
an atom in the dressed state \n) r is associated with a Berry connection 
A n (r) = ih(n\'V\n) r = ih (a* 1 Va n q + a* 2 Va n , 2 ) • This plays the role of 
a vector potential coupling to the motion of the atom. By careful choice of 
the spatial dependence of the fields, one can arrange that not only A n but 
also V X A n is nonzero (for example for the lowest energy dressed state), 
such that the particle experiences a nonzero artificial magnetic field [2i|. 
This method of generating artificial gauge fields has been demonstrated in 
experiments at NIST m using Raman coupling of spin states of rubidium, 
see the contribution by I. B. Spielman et al. 

For the method of Ref. [60] the vector potential is limited to \A\ < h/K 
where A r is the optical wavelength. By Stokes’ theorem, the total number 
of flux quanta through a circular atomic cloud of radius R is then = 
| / V x A d?r = j: <f A - dl < R/ A r . Hence, the flux density n^ = N^/^ttR 2 ) 
is limited to n<p < l/(i?A r ). In Ref. [61] it was shown that periodic optical 
dressing of internal states can allow the formation of “optical flux lattices”, 
with very much increased flux densities. These lattices lead to flux densities 
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ri 0 ~ 1 /Aj?, larger by a factor of R/ A r over previous proposals [2], which is 
R/ A r ~ 50 — 100 for typical experimental settings. 


Optical flux lattices 

The original motivation to consider periodic optical lattices coupling internal 
states [61] was to find a way in which to increase the strength of the artifi¬ 
cial magnetic fields to n$ ~ 1/A]:. Here we shall not focus on the real-space 
magnetic field. Instead, we concentrate on the band structure, defined in 
reciprocal space. We shall show that two-dimensional optical lattices involv¬ 
ing the coupling of internal states provide a very powerful way to generate 
topological (Chern) bands [62]. We shall return at the end to comment on 
the connection to real-space magnetic field. 

Consider an atom with N s long-lived internal states. We denote the (plane- 
wave) state of an atom of internal state a = 1... JV S and momentum hq 
by | a,q). The coupling of the atom to optical fields will lead to processes 
which change the internal state and/or the momentum of the atom. We 
denote the optical coupling for | a,q) —> \a',q') by V^'fl q . We consider the 
case of periodic lattices, for which the momentum transfers of the set of all 
such couplings, Vfl' a , are commensurate. The wave vector of any component 
a is then only conserved up to the addition of reciprocal lattice vectors 
G. Similarly, by Bloch’s theorem, the energy eigenstates can be assigned 
a band index n and a quasi-momentum hk , and decomposed as | i/> nk ) = 
E Q ,G C a k G\ a i k -9a-G) where g a accounts for possible momentum offsets 
of the different internal states 4 . The band energies E n (k ) follow from 


En( k K k G = 


9cc~ 


nk . \ ' taqq' 

GC aG + V 9ai _ 

cx',G ' 


+G'-g a -G C a k G' 


(1.24) 


where e q = h 2 \q\ 2 / 2 m is the kinetic energy for an atom of momentum hq. 

The structure of the couplings in Eqn (1.24) can be conveniently repre¬ 
sented by a lattice in reciprocal space. An example for N s = 3 is shown in 
Fig. 1.5 the sites denote the different values of ( a,G ); the bonds represent 
the couplings V£ a . (We provide here no experimental implementation of 
this lattice, but note that it can be generated by a variant of the triangular 
lattice in Ref. [ 63] .) 

It is instructive to construct the lowest energy band using this picture 
in the “nearly free particle” limit of small Consider first vanishing 


4 This labeling of Bloch states takes advantage of an enhanced spatial symmetry of the system 
und er co mbined translations and internal-state gauge-changes. For the model illustrated in 
Fig. |l.5^ a), the reciprocal lattice vector G\ is linked to a symmetry under the spatial 
translation by combined with the unitary transformation 

[|1)(1| +e i27r / 3 |2)(2| + e i47r / 3 |3)(3|]. 
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Figure 1.5 (a) Example of the reciprocal space representation of an optical 
lattice involving couplings of N s = 3 internal states. The sites denote the 
different values of (a, G ), and the bonds the couplings V^ a i+G ,_ g _ G . The 


offset momenta are g ± = 0 , g 2 = /ci, = 2 /ci, and the reciprocal lattice 
has basis vectors G\ = 3«i, G 2 = k 2 . (b) A representation of one plaquette 
of the reciprocal space lattice. For vanishing coupling V = 0, the lowest 
energy state is one of the plane-wave states ({|a), \b), |c), \d}}) defined in 
the text and is degenerate for quasi-momenta along the dashed lines. For 
nonzero V, gaps open along these lines. To first order in V, the inte gral o f 
the Berry curvature of the lowest band over this plaquette is by Eqn (1.26). 


coupling |V^ a | = 0. From Eqn (1.24), at a given quasi-momentum Hk, the 
lowest energy state is that state with smallest \k — g a — G\ 2 : that is, it 
is the plane wave state \a,k — g a — G) associated with the lattice point 
(a, G) closest to k. This divides each square plaquette into four regions, 
shown in Fig. |1.5[ b) for the leftmost plaquette of Fig. [E5| a), with the lowest 
energy state in these regions being |a) = |l,fc — gq), |6) = 12 , fc — g 2 ), 
|c) = \2, k — g 2 — G 2 ), | d) = 11, k — g 1 — G 2 ). The dashed lines in Fig. |1.5[ b) 
show the locations where the lowest energy state is degenerate for V = 0: 
two-fold degenerate on each dashed line and four-fold degenerate at the point 
k* where the dashed lines cross. 

For |V^ ,Q | 7 ^ 0, the degeneracies along the dashed lines are lifted. We 
describe the behaviour in the case where the couplings around this plaquette 
-Ve l<t> with = 4> a b,bc,cd,da as illustrated in Fig. 


are 


1.5 


b). Working to lowest 
order in V involves first order (degenerate) perturbation. At the point of 
four-fold degeneracy k* this leads to the Hamiltonian 


( 

V 


e k* 

_y 

0 

_y g i'&da 


_y g i&ab 

£k * 

_y g^^bc 


0 

_yg 'i&bc 


— y e i</>da \ 

0 

_yg i&cd 


0 


/ 


e k* 

y g^0cd 


(1.25) 
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The spectrum of this Hamiltonian depends on the gauge invariant phase 


= 4>ab + 4>bc + 4>cd + 4>da (1.26) 

For <hp 7 ! 7 r the lowest energy state is non-degenerate at k *, and indeed 
for all k within this plaquette. The lowest band therefore has a well-defined 
Berry curvature, with no singularities, for <f> p 7 ! it. One can then compute 
the integral of the Berry curvature over the plaquette in terms of the line 
integral of the Berry connection around the four sides of the plaquette. Still 
working to first order in V, the adiabatic transfer along each of these sides 
just involves the coupling of a pair of states (e.g. |a) and | 6 ) coupled by 
— Ve l ^ ab ). It is straightforward to show that the total integral of the Berry 
connection around the plaquette is just 4> p . Thus, to first order in V, the 
integral of the Berry curvature over the plaquette is precisely 4> p for $ p 7 ^ 7 r. 
For <h p = 7 T the lowest two bands touch at a Dirac point at k = k*, and the 
Berry curvature of the lowest band is ill-defined. 

This result provides a very powerful prescription by which to design op¬ 
tical couplings that generate Chern bands. By choosing the phases of the 
couplings V™' a around each plaquette in the reciprocal space lattice, one can 
specify (the integral of) the Berry curvature of the lowest energy band over 
each of these plaquettes. If the sum over all plaquettes in the first Brillouin 
zone <3?tot = SpeFBZ ^p * s nonzero, and all <h p / 7 r, then the lowest energy 
band will have a Chern number of C = ^> t ot/(27r). Thus the net number of 
“flux quanta” through this reciprocal space lattice sets the Chern number 
of the lowest energy band. For example, if the couplings of Fig. |1.5[ a) are 
chosen 


V 21 = V 32 = V 13 = -V 

v — K.1 — K-i ’ 


TyQQ t r i 27ra/3 

—«2 V e 


(1.27) 


then the lowest energy band will have integrated Berry curvature of 27 t/ 3 in 
each of the plaquettes in reciprocal space, so Chern number C = 1. Although 
this result applies in the weak lattice limit <C E r , the fact that the 

Chern number is a topological invariant guarantees robustness of the result 
up to moderate | Vjf “ |. 

Our considerations of the energy bands, and their topological properties, 
have been presented in reciprocal space. How do these considerations relate 
to the real-space magnetic field discussed in §1.4.2| ? The magnetic field is 
defined in the limit in which the particle moves adiabatically in real space , 
valid when |Vj? ,Q j E r . In this limit, the kinetic energy term in (1.24) 


is negligible and the energy eigenstates follow from a tight-binding model 
on the reciprocal space lattice. An elegant duality emerges [62] : the mag¬ 
netic flux through the unit cell in real space due to adiabatic motion in 
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the lowest energy dressed state is equal to the Chern number of the lowest 
energy band of the tight-binding model defined by the couplings V°' a on 
the reciprocal space lattice. To form an “optical flux lattice”, with non-zero 
magnetic flux for the lowest energy dressed state through the real-space 
unit cell, one should just choose the couplings V^' a on the reciprocal space 
lattice to give a lowest energy band with non-zero Chern number. For the 
example of Fig. |1.5[ a), with the above couplings (1.27), the reciprocal space 
model is precisely the Harper-Hofstadter model for flux 1/3 per plaquette. 
As discussed in Section 1.3.2| this model readily leads to Chern bands: for 
flux 1/3 the lowest band has unit Chern number, so the couplings (1.27) 
form an optical flux lattice with one magnetic flux quantum through each 
real-space unit cell. 

The non-zero flux densities of optical flux lattices lead to low-energy bands 
that are very similar to those of the lowest Landau level: with unit Chern 
number, and with narrow width in energy. The bands depart from the ex¬ 
act degeneracy of Landau levels owing to the fact that the magnetic held 
and the scalar potential are non-uniform in space. However, for well-chosen 
parameters these effects can be made small in practical implementations 


1.5 Discussion 

In this Chapter, we discussed how Chern bands can be created in cold-atom 
systems, through a wide variety of schemes based on atom-light coupling, ro¬ 
tation, or time-periodic modulations. Several of these techniques have been 
experimentally implemented in different laboratories, demonstrating physi¬ 
cal manifestations of synthetic magnetic fields M Em IM EH 03123 ESI E3 
fl5ll48j . Beyond the signatures discussed here, topological energy bands could 
also be probed in cold-atom systems through the detection of Skyrmion pat¬ 
terns in time-of-flight images [23, 25], interferometry |66l [67], or the presence 
of chiral edge-states [681 1MI [701 EH- 

This Chapter focused on single-particle phenomena in Chern bands. The¬ 
oretical studies have extended the ideas presented here in a variety of ways. 
The methods for generating Chern bands for cold atoms have been extended 
to higher dimensions and other symmetry classes, including the Z 2 topo¬ 
logical insulator in three dimensions [721113] both in lattice-based [THUS] 
and in continuum m formulations, and to systems with sublattice (chi¬ 
ral) symmetry [77] . It is very interesting to consider the effects of interac¬ 
tions on degenerate bosons or fermions in the flat, or nearly flat, topological 
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bands that we have described above. For the Harper-Hofstadter model, nu¬ 
merical studies have established the existence of fractional quantum Hall 
states, including the Laughlin state P3 E3 EH EU as well as states that 
exist only on the lattice [ 801 [82] . The use of internal states to form optical 
flux lattices leads to very flat topological bands which have been shown to 
give rise to exotic non-Abelian phases of bosons even for weak two-body 
repulsion [63 i 83] . as well as to interesting forms of ferromagnetic-nematic 
ordering for one-component fermions [EH] . Studies of fractional Chern insula¬ 
tors have identified a wealth of tight-binding lattice geometries with Chern 
bands where strongly correlated phases can appear mm and which may 
find realizations in cold atom set-ups [85] • These theoretical works show 
that cold atomic gases have the potential to provide an ideal setting in 
which to explore novel forms of strongly correlated quantum phases, includ¬ 
ing fractional Chern insulators and other lattice-based strongly correlated 
topological phases. 
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